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Abstract. An abstract framework for constructing stable decompositions of the spaces 
corresponding to general symmetric positive definite problems into "local" subspaces and 
a global "coarse" space is developed. Particular applications of this abstract framework 
include practically important problems in porous media applications such as: the scalar 
elliptic (pressure) equation and the stream function formulation of its mixed form, Stokes' 
and Brinkman's equations. The constant in the corresponding abstract energy estimate 
is shown to be robust with respect to mesh parameters as well as the contrast, which is 
defined as the ratio of high and low values of the conductivity (or permeability). The 
derived stable decomposition allows to construct additive overlapping Schwarz iterative 
methods with condition numbers uniformly bounded with respect to the contrast and 
mesh parameters. The coarse spaces are obtained by patching together the eigenfunctions 
corresponding to the smallest eigenvalues of certain local problems. A detailed analysis 
of the abstract setting is provided. The proposed decomposition builds on a method of 
Efendiev and Galvis [Multiscale Model. SimuL, 8 (2010), pp. 1461-1483] developed for 
second order scalar elliptic problems with high contrast. Applications to the finite element 
discretizations of the second order elliptic problem in Galerkin and mixed formulation, 
the Stokes equations, and Brinkman's problem are presented. A number of numerical 
experiments for these problems in two spatial dimensions are provided. 



1. Introduction 

Symmetric positive definite operators appear in the modeling of a variety of problems 
from environmental and engineering sciences, e.g. heat conduction in industrial (com- 
pound) media or fluid flow in natural subsurfaces. Two main challenges arising in the 
numerical solution of these problems are (1) the problem size due to spatial scales and (2) 
high-contrast due to large variations in the physical problem parameters. The latter e.g. 
concerns disparities in the thermal conductivities of the constituents of compound media 
or in the permeability field of porous media. These variations frequently range over several 
orders of magnitude leading to high condition numbers of the corresponding discrete sys- 
tems. Besides variations in physical parameters, the discretization parameters (e.g. mesh 
size) also lead to large condition numbers of the respective discrete systems. 

Since in general high condition numbers entail poor performances of iterative solvers 
such as conjugate gradients (CG), the design of preconditioners addressing the size of the 
problem and the magnitude in variations of the problem parameters has received lots of 
attention in the numerical analysis community. The commonly used approaches include 
domain decomposition methods (cf. e.g. [20 \ [271 [28] ) and multilevel/multigrid algorithms 
(cf. e.g. in HU [30] ) . For certain classes of problems, including the scalar elliptic case, these 
methods are successful in making the condition number of the preconditioned system 
independent of the size of the problem. However, the design of preconditioners that are 
robust with respect to variations in physical parameters is more challenging. 
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Improvements in the standard domain decomposition methods for the scalar elhptic 
equation, —V • {k{x)'V4>) = /, with a highly variable conductivity k{x) were made in 
the case of special arrangements of the highly conductive regions with respect to the 
coarse cells. The construction of preconditioners for these problems has been extensively 
studied in the last three decades (see e.g. |7l [HI [191 EOl USJ). In the context of domain 
decomposition methods, one can consider overlapping and nonoverlapping methods. It was 
shown that nonoverlapping domain decomposition methods converge independently of the 
contrast (e.g. \18\ |2T] and [28i Sections 6.4.4 and 10.2.4]) when conductivity variations 
within coarse regions are bounded. The condition number bound for the preconditioned 
linear system using a two level overlapping domain decomposition method scales with the 
contrast, defined as 

maxK(a;)/minK(a;) (1-1) 

where $7 is the domain. The overall condition number estimate also involves the ratio H/6, 
where H is the coarse-mesh size and 5 is the size of the overlap region. The estimate with 
respect to the ratio H/5 can be improved with the help of the small overlap trick (e.g. 
|14l [28] ) . In this paper we focus on improving the contrast-dependent portion (given by 



(1.1)) of the condition number. A generous overlap will be used in our methods. 

Classical arguments to estimate the preconditioned condition number of a two level 
overlapping domain decomposition method use weighted Poincare inequalities of the form 



/ K\V^\'^{ij-I^'iljfdx<C I K\V^l)\'^dx, (1.2) 



where is a local subdomain in the global domain 17, ^ is a partition of unity function sub- 
ordinate to Li;, and '\\) E H^{u}). The operator I^ip is a local representation of the function 
■0 in the coarse space. The constant C above appears in the final bound for the condition 
number of the operator. Many of the classical arguments that analyze overlapping domain 
decomposition methods for high contrast problems assume that '^i'"} is bounded. 



When only |V^| is bounded, inequality (1.2) can be obtained from a weighted Poincare 



inequality whose constant is independent of the contrast. This weighted Poincare inequal- 
ity is not always valid, so a number of works were successful in addressing the question 
of when it holds. In [26], it was proven that the weighted Poincare inequality holds for 
quasi-monotonic coefficients. The author obtains robust preconditioners for the case of 
quasi- monotonic coefficients. We note that recently, the concept of quasi-monotonic co- 
efficient has been generalized in |22| [23] where the authors analyze nonoverlapping FETI 
methods. Other approaches use special partitions of unity such that the "pointwise en- 



ergy" K(ic)|V^p in (1.2) is bounded. In [Ill[29], the analysis in [28] has been extended 
to obtain explicit bounds involving the quantity k{x)\\/^\'^. It has been shown that if the 
coarse space basis functions are constructed properly, then «;(a;)|V^p remains bounded for 
all basis and partition of unity functions used, and then, the classical Poincare inequal- 
ity can be applied in the analysis. In particular, two main sets of coarse basis functions 
have been used: (1) multiscale finite element functions with various boundary conditions 
(see [U [T^ 117] ) and (2) energy minimizing or trace minimizing functions (see \29\ [33] 
and references therein). Thus, robust overlapping domain decomposition methods can be 
constructed for the case when the high-conductivity regions are isolated islands. These 
methods use one coarse basis function per coarse node. 

Recently, new coarse basis functions have been proposed in |10l [TT] . The construction 
of these coarse basis functions uses local generalized eigenvalue problems. The resulting 
methods can handle a general class of heterogeneous coefficients with high contrast. A 
main step in this construction is to identify those initial multiscale basis functions that 
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are used to compute a weight function for the eigenvalue problem. These initial multiscale 
basis functions are designed to capture the effects that can be localized within coarse 
blocks. They are further complemented using generalized eigenfunctions that account 
for features of the solution that cannot be localized. The idea of using local and global 
eigenvectors to construct coarse spaces within two-level and multi-level techniques has 
been used before (e.g. [HI [25]). However, these authors did not study the convergence with 
respect to physical parameters, such as high contrast and physical parameter variation, and 
did not use generalized eigenvalue problems to achieve small dimensional coarse spaces. 

In many applications, the discretization technique is chosen so that it preserves the 
essential physical properties. For example, mixed finite element methods are often used 
in flow equations to obtain locally mass conservative velocity fields. In a number of flow 
applications, high-conductivity regions need to be represented with flow equations, such 
as Navier-Stokes' or Stokes' equations due to high porosity. Such complex systems can be 
described by Brinkman's equations that may require a special stable discretization. The 
method proposed in TIJ cannot be easily applied to vector problems and more com- 
plicated discretization methods. More precisely, it requires a proper eigenvalue problem 
for each particular differential equation. 

In this paper, we extend the framework proposed in |10| [TT] to general symmetric 
bilinear forms. The resulting analysis can be applied to a wide variety of differential 
equations that are important in practice. A key in designing robust preconditioners is 
a stable decomposition of the global function space into local and coarse subspaces (see 
|28j). This is the main focus of our paper. We develop an abstract framework that allows 
deriving a generalized eigenvalue problem for the construction of the coarse spaces and 
stable decompositions. In particular, some explicit bounds are obtained for the stability 
constant of this decomposition. In the scalar elliptic case, the analysis presented here leads 
to generalized eigenvalue problems that differ from those studied in jlOl [TT] . 

We consider the application of this abstract framework to Darcy's equation as well as to 
Brinkman's equations. Brinkman's equations can be viewed as a generalization of Darcy's 
equation that allows both Darcy and Stokes regions in the flow. Because Darcy's equation 
is obtained under the assumption of slow flow, Brinkman's equation is inherently high- 
contrast. It combines high flow described by Stokes' equations and slow flow described 
by Darcy's equation. In fact, the use of high conductivities in Darcy's equation may be 
associated to free flow or high porosity regions that are often described by Brinkman's 
equations (see [5]). The proposed general framework can be applied to the construction 
of robust coarse spaces for Brinkman's equations. These coarse spaces are constructed by 
passing to the stream function formulation. 

In the paper at hand, we also discuss some coarse space dimension reduction techniques 
within our abstract framework. The dimension reduction is achieved by choosing initial 
multiscale basis functions that capture as much subgrid information as possible. We discuss 
how the choice of various initial multiscale basis functions affects the condition number of 
the preconditioned system. 

To test the developed theory we consider several numerical examples. In our first set 
of numerical examples, we study elliptic equations with highly variable coefficients. We 
use both piecewise bilinear and multiscale basis functions as an initial coarse space. In 
the latter choice, they capture the effects of isolated high-conductivity regions. Our nu- 
merical results show that in both cases one obtains robust preconditioners. The use of 
multiscale basis functions as an initial coarse space allows substantial dimension reduction 
for problems with many small isolated inclusions. Both Galerkin and mixed formulations 
are studied in this paper. The next set of numerical results are on Brinkman's equa- 
tions. These equations are discretized using /7(diw)-conforming Discontinuous Galerkin 
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Figure 2.1. Triangulation 7^. 



methods. The numerical results show that the number of iterations is independent of the 
contrast. In our final numerical example, we consider a complex geometry with highly- 
variable coefficients without apparent separation of high and low conductivity regions. 
The numerical results show that using multiscale initial basis functions, one can obtain 
robust preconditioners with small coarse dimensional spaces. 

The paper is organized as follows. In Section [2j problem setting and notation are intro- 
duced. In Section |3} the abstract analysis of the stable decomposition is presented. Sec- 
tion |4] is dedicated to applications of the abstract framework to the Galerkin and mixed 
formulation of Darcy's equation. Stokes' equations, and Brinkman's equations. In Sec- 
tion [5| we discuss the dimension reduction of the coarse space. Representative numerical 
results are presented in Section [6j 

2. Problem Setting and Notation 

Let C be a bounded polyhedral domain, and let Th be a quasiuniform quadrilateral 
(n = 2) or hexahedral {n = 3) triangulation of il. with mesh-parameter H. Let X = 
{^j}j=i be the set of nodes of Th, and for each Xj G X we set 



interior 



[J{T\TgTh, x,gT} 



i.e., i}j is the union of all cells surrounding xj (see Figure 2.1). 

For a suitable separable Hilbert space 'fo = i^o^^l) of functions defined on Q and for any 
subdomain u C Q we set '^(w) := {^Itjl^ G '^o} and we consider a family of symmetric 
positive semi-definite bounded bilinear forms 

For the case a; = (7 we drop the subindex, i.e., a (•, •) = oq (•, •) and we additionally 
assume that a (•, •) is positive definite. For ease of notation we write a^j {4>, tp) instead of 
o-u) {4'\u}, V'U) for any (/>, -0 £ ^o- Furthermore, we assume that for any (j) £ Yq and any 

family of pairwisely disjoint subdomains {uj}^^-^ with = ^ 



« = X] 

Now, the goal is to construct a "coarse" subspace 
property: for any cj) £ i^Q there is a representation 



(2.1) 

Tni^) of % with the fohowing 



rial 

j=0 



with 00 £ i^H, and cpj £ ^^{flj) for j = 1, 



(2.2a) 
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such that 

^a((/>,-,(/.,)<Ca (</),</.), (2.2b) 

j=0 

where '^(ilj) C is a subspace such that a^^ (•, •) : {'^{^Ij), '^(ilj)) — )■ M is positive 

definite. Note that uq. (•, •) : {Y(Qj),'f{0,j)) — t- M is in general only positive semi- 
definite. 

By e.g. j28l Section 2.3] we know that the additive Schwarz preconditioner corresponding 



to (2.2a) yields a condition number that only depends on the constant C in (2.2b) and the 
maximal number of overlaps of the subdomains j = 1, . . . , rix- Thus, we would like to 
"control" the constant C and keep the dimension of "as small as possible". 

For the construction of i^n we need more notation. Let {^j}"^]^ : — )• [0, 1] be a 

partition of unity subordinate to {Qj}^^-^ such that supp(^j) = ilj and for any (p £ 
we have ^jcj) G ^o, {(,j4')\nj G j = ^, ■ ■ ■ ,nx. Using this notation we may for any 

i, j = 1, . . . ,nx define the following symmetric bilinear form: 

mn,{4>,'4') ■■= "^an^ {^j^i4>, ^j^iip) ^'^'^^ 
where Ij := {i = 1, . . . ,nx\^i D ^ 0}. Let rij := max ^Ij denote the maximal 

j = l,...,?la, 

number of overlaps of the fij's. To ease the notation, as we did for the bilinear form 
a (•, •), we write mn^ ((/>, ip) instead of mn. {(pluj, V'ln,) for any (p, ip G %. 



Due to our assumptions on we have that (2.3) is well-defined. Also note, that 



since supp(^j) = J7j we have £,j4>\nj =0 4^ (f)\n. = 0, which implies that mQ. (•, •) : 
M is positive definite. 

Now for j = 1, . . . , ?7-a; we consider the generalized eigenvalue problems: Find (A^, (p^) G 
(M, y{nj)) such that 

an, (V, ^i) = A>n, {ip, ^i) , G (2.4) 

Without loss of generality we assume that the eigenvalues are ordered, i.e., < A{ < A2 < 
...<A|<Al+,<.... 

It is easy to see that any two eigenfunctions corresponding to two distinct eigenvalues are 
o-Uj (■) ■)" ^iid niQj (•, •)-orthogonal. By orthogonalizing the eigenfunctions corresponding 
to the same eigenvalues we can thus, without loss of generality, assume that all computed 
eigenfunctions are pairwisely an- (•, •)- and (•, •)-orthogonal. Now, every function in 
"/(Qj) has an expansion with respect to the eigenfunctions of ( |2.4| ). This is the reason 
why the generalized eigenproblem is posed with respect to as opposed to Yo^Qj). 

For (p G let (pQ be the m^^. (•, •)-orthogonal projection of <P\q, onto the first Lj 



eigenfunctions of (2.4), where Lj G No is some non-negative integer, i.e., 

mnJ(p-cplip{)=0, yi = l,...,Lj. (2.5) 



If Lj = 0, we set <pl = 0. We assume that for any j = 1, . . . ,nx and a sufficiently small 
"threshold" > we may choose Lj such that ^^.^i > '''x^- Since any function in 



Y{Qj) has an expansion with respect to the eigenfunctions of (2.4) we can now choose 
a, S M so that 
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Then we observe that 



\i>Lj i>Lj J 

= ^ m^^ ("jW' '^i'^i) (by orthogonahty) 



i>Li 



^an, [ai^i, (by ^ 



(2.6) 



i>Lj i 



< Tx an, 



3. Coarse Spaces yielding Robust Stable Decompositions 



With these prehminaries we are now able to define a decomposition described in (2.2): 
First, we specify the coarse space by 

i^H ■■= span{^jipl I j = 1, . . . , Tia, and i = 1, . . . , Lj}. (3.1) 

Then, for any (p € Y let 

ria! 

h-=Y.^3<l^o^^H, (3.2a) 



where (f)^ is chosen according to (2.5). For j = 1, . . . ,nx define 

^3 ■■= i^ji^ - <Po))\n, e >o(%) so that <j) = Y<Pj- (3.2b) 

j=0 

Before analyzing this decomposition we summarize all assumptions using the notation 
above: 

(Al) a(j (•, •) : ('^((x>), 'V{uj)) — 7- M is symmetric positive semi-definite for any subdomain 

w C il, and a (•, •) = a^ (•, •) is positive definite. 
(A2) For any (p £ and any pairwisely disjoint family of subdomains {ujjYj^i with 

I^LJj = we have a (</>, <f)) = Y^'^'^i ('^' (This implies that a (0, 0) < 

i=i i=i 

^^an, (i?^, 1?^*), since (•, •) is semi-definite.) 

i=i 

(A3) an,. (•, •) : M is positive definite for all 1 < j < n^,. 

(A4) {i^j}"^^ : — )• [0, 1] is a family of functions such that 

(a) X] = ^ ^' 
i=i _ 

(b) supp(^j) = VLj for j = 1, . . . , na;; 

(c) For G ro we have ijCp G and (Cj(A)b, G ■^o(^) for j = 1, . . . , n^. 

(A5) For a sufficiently small threshold r^^ we may choose Lj such that A]^ ,_|_j^ > t^^ for 
ah j = 1, ... ,713,. 
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As noted above, these assumptions imply in particular, that for any 1 < j < Ux, the 
bilinear form m^^. (•, •) : — )• M is positive definite. 

The following lemma establishes the required energy bound for the local contributions 
in decomposition (|2.2|): 



(3.3) 



Lemma 3.1. Assume (Al)~(A5) hold. Then, for (j) & 'Vq we have 



i=i 

where C only depends on nj (the maximal number of overlaps of the subdomains {Qj}^!^) 
Proof. Observe that 



(by (3.2b)) 



(by (3.2a)) 



i=l 



i=l 



^ V ' 

=.Ei 

where in the last step we have used Schwarz' inequality together with fflj < nj. Note, 
that we furthermore have 

^1 = Y.Y.''^d^M<p-^h),m 



j=i ieij 



i=l jeh 

Tlx 

^mn, {<P-4>1, (t}-(t>o) 



(by (2.3)). 



i=l 



Thus, we obtain 



Tlx '''X 

E " ('^i' "^j) - '^'fE '^^j -(t^o^4>- 0o) 



< n E 

< n]Txa ((/), (j)) , 



(by ^ 



(3.4) 



where the last inequality holds due to (A2) and the fact that any point in Vt simultaneously 
lies in at most rij subdomains J7j , j = 1 , . . . , ria; . □ 



Remark 3.2. Note, that by the proof of Lemma 3.1 we in particular have 
Tlx I \ 

E ^^E ^^('^ - '^o), e^E ^^('^ - '^o) < n] TA a (</), <t^) . 
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Now, we proceed with the necessary energy estimate for the coarse component <j)Q of the 



decomposition (3.2a). 



Lemma 3.3. Let (A1)-(A5) hold. Then, for (po defined by (3.2a) we have that 

a (00, M < (2 + CrA)a(0, (j)) , 
where, as above, C only depends on nj. 
Proof. First, we note that 



(3.5) 



\i=l 1=1 
i=l i= 

2a 5^^,(0^ -./.), 



< 



vi=l i=\ 



.(00 



+ 



+2a 



(by ^2^ ) 
(by (A4)) 



(3.6) 



--■■E2 



where we have used Schwarz' inequahty. Now, observe that 
E2 < Y^an^ [^e.(<^^_<^), 



(by (A2)) 



I (1 - 6- + Ci)E ^^('^0 - 0), (1 - + ^,)E ^^(-^0 



< 2^ an, (l-e,)E^^('^o-0), (l-e.)E^«('^o 

j=i \ iei, leij 

^ V 

=:E3 

rial I 

< 2^3 + 2n2 TA a (0, 0) , 
where we have again used Schwarz' inequahty. Note, that 



^3 = E«^. E ^'^^^('^o 

j=i Vze/Ab} ie-f. 



< n,E E U;E^*('^o-0),e/E^«('^o 

j=i ieij\{j} 



(by Remark 3.2) 
(3.7) 



(by (A4)) 
(by Schwarz' inequahty). 



Since 



an, |e/E^'('^o-0),aE^^('^o 

is/,- ie/,- 



= «f^i 6 E ^»('^o-0),6 E ^^("^0 

< aj,, (6E^^('^o-0),ezE^*('^o-' 
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we thus have 



^3 < njY, E ^nA^lY.^^{<l^i-c|>),^lY,u^l> 

rise I 
rial I 

1=1 \ ieii ieii 
< njrxa {(p, </>) 



(by Remark [3^. 



Combining (3.6), (3.7), and (3.8) we obtain 

a (00, 0o) < (2 + 4{nj + n'j)Tx) a 

Combining Lemmas |3.1| and |3.3| we have the following: 



(3.8) 



□ 



Theorem 3.4. Assume (A1)-(A5) hold. Then, the decomposition defined in (3.2) satisfies 

rial 

Y,a{^j,^j)<{2 + CTx)a{<P,4>), (3.9) 

3=0 

where C only depends on nj. 

4. Applications 

To apply the theory developed in Sections [2] and [3] to some particular problem we need 
to verify assumptions (Al)-(A5). Once this is established, we can conclude that the 
condition number of the corresponding additive Schwarz preconditioned system has an 
upper bound which depends only on n/ and t\. Thus, if we can show the validity of (A5) 
with T\ chosen independently of certain problem parameters and mesh parameter the 
condition number will also be bounded independently of these problem parameters and 
the mesh parameter. 

In addition to verifying assumptions (Al)~(A5) we also have to make sure that the 



number of "small" eigenvalues in (2.4), i.e., those below the threshold r, , is a small in 



order for our method to be practically beneficial. Note, that the choice of ^ > and thus 
the number of small eigenvalues is not unique. Nevertheless, for a certain choice of r^^ 
and a given problem we may still aim to establish that the number of eigenvalues below 
the chosen threshold is uniformly bounded with respect to changes in specific problem 
parameters. In this case the additive Schwarz preconditioner corresponding to the stable 



decomposition (3.2 ) has a coarse space whose dimension is uniformly bounded with respect 



to these parameters. 

4.1. The Scalar Elliptic Case — Galerkin Formulation. As a first application of the 
abstract framework developed above, we consider the scalar elliptic equation 

- V • (kVc/)) = /, :e G il, and = 0, dVt, (4.1) 

where k G L°°(r2) is a positive function, which may have large variation, G /7g(i7), 
and / G L^(r2). Note that we can always reduce to the case of homogeneous boundary 
conditions by introducing an appropriate right hand side. 
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Figure 4.1. Subdomain with connected components of 0^. In the present 
configuration Lj = 7. 



With 'fo := Hq{Q,), the variational formulation corresponding to (4.1) is: Find (j) £ 'fo 
such that for all ip €z i^o 

a^^{(j),i}) := f K{x)V(t)-Vi;dx= [ f^jdx. 
Jn Jn 

Let Yoi^j) '■= -f^o (^j) C Yolnj for any j = 1, . . . ,nx. Choosing the Lagrange finite 

element functions of degree one corresponding to Th, we readily see that for the scalar 
elliptic case, i.e., when setting a (•, •) = a^^ (•, •), assumptions (A1)-(A4) are satisfied. 
Let us for now assume that k assumes only two values. More precisely, 

k{x) = < ^ with $7* U = $7 and Kmax ^ '^min > 0. 

[ Kmax in OP, 

Without loss of generality, we may take Kmin = 1- Now we wish to establish the existence 
of such that the number of eigenvalues below this threshold is independent of the 
contrast Kmax/Kmiri! which is the problem parameter of interest in this situation. 



By the well-known min-max principle [21], we know that the i-th. eigenvalue of (2.4) is 
given by 

^3 ■ «nf(V',V') 

A;; = mm max — -, (4.2) 

where 'fi{Vtj) is any i-dimensional subspace of (Here mf^ (•, •) is defined according 

to with a (•, •) replaced by a^^ (•, •)•) 



We denote fi^ := il^nOj and := QPr\Q.j. Now, for each i = 1, . . . , let Q| ^ denote 
the i-th connected component of ilj, where Lj denotes the total number of connected 
components of O^. If 0^ = 0, we set ^ = ilj and Lj = 1. 
Now, we define the following subspace of '^{Vtj): 

ysEi^j) ■= I*/- e ^i^j) I y ^ (/'rfa; = , for i = 1, . . . . (4.3) 
It is straightforward to see that any Lj • + 1-dimensional subspace of 'f'[Vtj) has a non-trivial 



intersection with '^g^i^j)- Thus, by (4.2) there exists a non-zero (p S '^sE^^i) such that 



alf ((/), ct>) 
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Using the definitions of (•, •) and m.^ (•, •) we note that 



ielj ^1 

< 2^ / H(v(eiei)'/'i' + ioe*v0p) 

< Cy [ K((H-^^f + \V(l)\^) dx 



(4.5) 




< Cni \ H-'^ / K(t)^dx + af^^ 




where C only depends on the choice of the partition of unity. Furthermore, we observe 
that 

i Kcj)^ dx = I (jP' dx + KmaxT^] / 0^ dx 

< CH^ ij jVc^\^dx + K^^J2 J N^^l^dx 

\ i=l ^Jii J 

= CH"^ 1 K\V(t)\^dx, 

where we have used Poincare's inequality, which is possible since (p S '^sEi^j)- Here, C is 
a constant which only depends on the geometries of and ^ji, i = 1, . . . ,Lj. Thus, we 
obtain 

m^^ (</., 0)<Ca^^(<A, 0), 



which together with (4.4) yields a uniform (with respect to Kmax/'^min and H) lower bound 
for ^L.^i- Thus, we have verified assumption (A5) with tx independent of Kmax/^min and 
H. In particular we see that for a suitably chosen t\ the number of generalized eigen- 
values below and satisfying (2.4) is bounded from above by the number of connected 
components in il.'j, i.e., Lj. 

4.2. The Scalar Elliptic Case — Mixed Formulation. In this Subsection we consider 
the mixed formulation of the scalar elliptic equation, also known as Darcy's equations, in 
2 spatial dimensions, i.e., n = 2, modeling flow in porous media 

= / in ri, 

= inn, (4.6) 

= on on. 

Here, p G Lo(^) •= L'^i^)/^ denotes the pressure, u G H{div; fl) := {v € (^^(rj))^ | V • 
V G L^(i7)} is the velocity, / S (L^(r2))^ is a forcing term, and n denotes the unit 
outer normal vector to dO,. The viscosity is a positive constant, and k G L°°{Q) is a 
positive function. With Ho{div; 0,) := {v € H(div; i})\v ■ n = on dO,} the variational 
formulation of Darcy's problem is given by: Find {u, p) G {Ho{div; i7), Lq{Q)) such that 
for ah {v, q) e {Ho{div; Q), -^^o(^)) h.a,Ye 

/ fiK^^u ■ V dx — / pV ■ V dx — / qV ■ udx = / f ■ v dx. (4.7) 
Jn Jn Jn Jn 
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It is well-known (see e.g. |4j Chapter 12, p. 300]) that problem (4.7) is equivalent to the 
following problem: find u in the subspace Ho{divQ; Q) := {v G HQ^div] ^2) | V • = 0} 
such that 

/ ^k'^u ■ vdx = f ■ V dx G Hf)[divo; 
Jn Jn 

Let us additionally assume that 17 is simply connected. Then, we have (see e.g. |12j ) 
that there exists exactly one S Hq{Q) such that Vxcf) = u, where Vxcf) := — . 

This leads to the variational form of Darcy's problem in stream function formulation: Find 
(f) ^ := Hq{Q) such that for all ■0 G we have 



, -0) := / fiK ^Vx(j) • Vxip dx = / f-{Vx-ip)dx. 



(4.8) 



Note that a {(f), tp) = / /^k Vi;^ • Viljdx. Thus, for '^(fij) and as chosen in Subsec- 
tion 4.1 for j = 1 



,nx we can readily verify assumptions (A1)-(A4). Let us assume 
that, as above, k only assumes two values. Then, we can perform exactly the same ar- 
gument as in Subsection 4.1 (with k replaced by f^K~^) to establish the validity of (A5). 
This in turn establishes the robustness (with respect to K^a.x/ t^mm and H) of the stable 



decomposition (3.2) corresponding to the stream function formulation with its bilinear 
form (•, •). Thus, we may robustly precondition (|4.8|), namely, solve for (f> and recover 
u 



Vx(j) from (4.7) 



An equivalent approach, that we use in Section |6| to compute a solution of (4.7), is 
somewhat different and outlined in the following Remark. 

Remark 4.1. Instead of solving the stream function formulation for (j) and then recovering 
'Vx(p, one may equivalently use the coarse space corresponding to (4.8) to construct 



u 



a coarse space corresponding to (4.7) by applying Vx to the coarse stream basis functions. 



This then yields an equivalent robust additive Schwarz preconditioner for (4.7) (for details 
see ^ Section 10.4.2]). 

4.3. Stokes' Equation. As for the mixed form of the elliptic equation we assume that 
C is simply connected. Then we consider Stokes' equations modeling slow viscous 
fiows 



-/uAu + Vp 
V • u 



u 



f 





in 0, 
in 17, 
on 317, 



where p G LI{Q), u G (i?o(r2))^, / G (^^(17))^, and /u G M+. The variational formulation 
of the Stokes problem is: Find (w, p) G {{Hl{Q))'^ , -^o(^)) such that for all {v, q) G 
{{Hl{n)f, L2(17)) we have 



/iVm : Vv dx — I pV ■ v dx — / q\7 ■ u dx 



f ■ V dx, 



(4.9) 



where Vii : Vir := t^^tt^ denotes the usual Frobenius product. 



dx j dx 



Analogously to Section [4^ we can formulate an equivalent problem for stream functions: 
Find G ^0 :={?/; G ^2(17) n ^0^(17) | ^{qq = o} such that for all (j) e 

a'^ {4>, ip) := fiV{Vx(f)) ■.V{Vx^)dx= f-Vxtpdx. (4.10) 
Jn Jn 
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For a sufficiently regular partition of unity {Ci}"=i and spaces 'fo{il.j), j = l,...,nx, 
defined as %{nj) := jV' G n HliQj) \ |^|af^ . = o}, we can readily verify (Al)- 

(A4). For the verification of (A5), we define for j = 1, . . . ,nx 



cl)€r{^}j)\ (t) = 0, / V(l)dx = 



As above, it is straightforward to see that any 4-dimensional subspace has a non-empty 
intersection with Thus, by again using the min-max principle we see that there 

exists a (p £ '^si^j) such that 



rrir 



(4.11) 



Using the definitions of Uq. (•, •) and m^, (•, •) (see (2.3) with a (•, •) replaced by a (•, •)) 
we note that 



< Cni [ fi{4>'^H-^ + n{(p) ■.n{^) + H-^{V(t>f)dx, 
where C only depends on the particular choice of the partition of unity. Here T-L{4>) 

d^6 



dxidxj 



denotes the Hessian of 



the tensor product, and 



denotes the 



Frobenius norm associated with the Frobenius product defined above. Since (p G "V^{Vtj) 
we may apply Poincare's inequality to (p and its first derivatives. Thus, we obtain 



. (</., 4>) < Cm [ ii-H{(p) : n{(P) dx = Ca^^ (</., 0) , 
Jn, 



(4.12) 



where C only depends on the partition of unity and the shape of Qj but is independent 
of H. Combining (4.12) with (4.11) we obtain A4 > C with C independent of H. This 



verifies (A5) and thus establishes that the decomposition (3.2) and the corresponding 



additive Schwarz preconditioner are robust with respect to H. As outlined in Remark 4.1 



we can also obtain a robust additive Schwarz preconditioner for (4.9). 



4.4. Brinkman's Equation. As for Darcy's and Stokes' problem, we assume that ^2 C 
is simply connected. Brinkman's problem modeling fiows in highly porous media is given 
by (cf. [5]) 

—nAu + Vp + fiK^^u = f inQ, 



V • u 
u 



in 0, 
on 



(4.13) 



where p, /, and /i are chosen as in the Stokes' case and k as in the Darcy case. The 
variational formulation of the Brinkman problem is: Find {u, p) £ ((i?Q (il))^, Lq(J7)) such 
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that for all {v, q) G {{Hl{VL)f, LI{VL)) we have 



jiVu : Vv dx + 



fiK ■ V dx 



p\7 ■ V dx — / gV • u dx 



Jn 



dx. (4.14) 



Again, we adopt the setting of stream functions. For 'Vq as in section (4.3), the vari- 
ational stream function formulation reads: Find 
have 



G ii) such that for all -0 ^ % we 

a^((/>, V):= / Ai(V(Vx(/.) : V(Vx0) + K-^Vx0- VxV') da; = f-Vxt/jdx. (4.15) 
Jn Jn 

With and '^(ilj) as in section 4.3 for j = 1, . . . , n^, we readily verify (A1)-(A4). 
Note that 

((/,, t/j) = (0, -0) + (</>, Tp) and (</>, -0) = {4>, V') + m'^ (</>, "0) > 

where (•, •) is defined according to (2.3) with a (•, •) replaced by (•, •). 

Since for any if) G '^{VLj), j = 1, . . . ,nx, we have m^. (-;/', V') > (V') V') > and 



a^, (-0, 0^) , (V'; 0^) > we have 



(V-, 0^) 



<. (V', V') + <. (V-, 0) 

> mm 



(V', V') 



m^ . {ip, ip) m^, (?/;, V) 



(4.16) 



m^^ (ip, i)) mf^^ (V', V) + mg^. V-) 

This is an immediate consequence of the following inequality, valid for /3i , > and 

/33, /34 > 



^1 



/3i+/32 

/33+/34 



+ 



/33/34 



/33/34 

1 1 



> < 



Pi 

I /33^ 



' /33 - /34 

if^>^ 

/34 - /33 



> min 



/3i /32 

/33' /34 



4.2 



Combining (4.16) with the results from Subsections 4.2 and 4.3 we obtain 

■^max{Lj+l, 4} - ^' 

where C is independent of H and Kmax/'^mim and where Lj is chosen as in Subsection 
This verifies (A5) and thus establishes that the decomposition introduced by (3.2) and 
the corresponding additive Schwarz preconditioner are robust with respect to H and 
'^max/'^mm- As for the Darcy and the Stokes case, we can also obtain an equivalent robust 
additive Schwarz preconditioner for (4.14) (see Remark 4.1). 



5. Reducing the Dimension of the Coarse Space 

In the exposition above, no assumptions except (A4) and (implicitly) (A5) were made 
about the choice of the partition of unity In this section, we investigate possi- 

bilities of making a particular choice of the partition of unity denoted by {^j}"^^ that 
results in a reduction of the dimension of the coarse space Yh- The idea is that by re- 
placing {0}"=! with {Cj}^^ ^ one can avoid the asymptotically small eigenvalues for those 
connected components of Q.'^ which do not touch the boundary of any coarse cell T G Th- 
For this, we again assume the scalar elliptic setting as in Section [4. 1[ 

Let be a standard partition of unity as above. For each j = 1, . . . , ria, and each 

r C r^j, let be a solution of 

-V • (kV^j) = 0, in T, Ij = ij, on dT. 
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In- □ 



rijj^, k — I, . . . , Lj 



ni,, k = Lj + i 



, . . . , Lj 



Figure 5.1. Subdomain with connected components of ilj. The black 
connected components are those touching an edge interior to Qj. The grey 
connected components are the remaining ones. In this configuration Lj = 4 



and Lj = 7. 



The corresponding variational formulation reads: Find ^jlr G Hq{T) such that for all 
^p G Hq(T) we have 

4''fe,V')=0, yijeH^iT). (5.1) 



We set = in and check whether with {Cj}"=i instead of (A4) and (A5) 

are satisfied. 

As in Subsection 



4.1 



let , k 



1, . . . ,Lj be the fc-th connected component of Qj. 
Without loss of generality we may assume a numbering such that ^ for k = 1, . . . , Lj{< 
Lj) is a connected component for which ^'j {dT\dO,j) ^ for some T G il.j. Note, that 
in general Lj < Lj and that Lj — Lj is precisely the number of connected components of 
Qj which do not touch an edge of T C interior to Qj, i.e., for i = Lj + 1, . . . , Lj we 

have that 17^^ n {dT\dUj) = 0. We also define Qj := ^j\{ These notations 

k=Lj+l 



are illustrated in Figure 5.1 



To proceed we make the following assumption: 
(A) For S,j defined as above, we have 

< CH^^ and 



where C is independent of Kmax/'^min and H. 



Note that by [9„ Lemma 3.1] 



< CH" and that by essentially the same 



argument as in 1141 Theorem 4.3 and 4.51 we have that KmaxVf,- _ < CH'^ ^, 

^ * J max Sj LHn,\n,) " 

where as above n denotes the spatial dimension. To the best of our knowledge, obtaining 
rigorous L°°-estimates as stated in (A) is still an open problem and beyond the scope 
of this paper. Nevertheless, we consider reducing the verification of (A4) and (A5) to 
establishing (A) a significant improvement. 

It is easy to see that = 1 and that supp(^j) = flj for any j = 1, . . . ,nx. Thus, to 

establish the validity of (A4) it remains to verify that (^jip £ = HQ{fl) and {(,j'ip)\n- £ 
'f(){^lj) for all ip G ii). For this we restrict to the case of two spatial dimensions, i.e., n = 2: 
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Note that for some e > we have that £ H^~^''{Qj) (cf. ^1^)- Thus, by [1^ Theorem 
7.57] we know that S L°°(i}j). Using this and (A) we see that 



L2(nj) 



< 



+ 



H 



< oo. 



It is furthermore easy to see that supp(^j'i/') C Jlj, which estabhshes (A4). 
Similarly to section 4.1 we now define 



^SEi^j) ■■= {(t>^ I / (l)dx = Q ,ioi k = l,...,Lj 



(5.2) 



and by the min-max principle we know that there exists a </> G ^sEi^j) such that 



Ai > 



4'' {<t>, 



\SE , 



(5.3) 



where mf^ (•, •) is defined as mf^ (•, •) with {^j}"^;^ replaced by {^j}^^]^. 

In order to obtain a uniform (with respect to Kmax/ ^min and H) lower bound for A~ 

Lj + l 

and thus establish (A5), we need to verify that 

m^f((/<,0)<Ca^f ((/>,</.), (5.4) 
with C independent of Kmax/'^min and H. 



Revisiting estimate (4.5), we obtain 



mgf {ct>, 0) < 2^ / Aj|V(Cj6.)0l' dx + 2a 



SE 



Thus, it suffices to bound for any i G 



K\v{ijiiWdx= I K\v{ijii)(i,\'dx+ /^^__ /t|v(e,-e*)<Ardcc 



■:Ei 



by af^ {4>, (/)). To avoid unnecessary technicalities, we make the simplifying assumption 
that each connected component of flj contains at least one ilj ^ with k = 1, . . . ,Lj. If 



this assumption is violated, one simply needs to introduce additional conditions in (5.2) 
ensuring that the average of functions is zero over each connected component of that 
does not contain any ^ with k = 1, . . . ,Lj. 

Assuming (A), we can find the required estimate of E/j^ as follows: Proceeding as in 
Subsection |4JJ we see 



^4 < C H-'^ [ K(t>^dx<C [ k| V0|2 dx<C a, 



SE , 



where we have used Poincare's inequality, which is possible since (p G "^sE^^i)- to (A), 
C can be chosen independently of Kmax/ /tmin and H, but it may depend on the geometry 



of ilj. 
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For an estimate of note that by (A) 



< 



, 



where we have used Poincare's inequahty. This estabhshes (5.4), which yields the vahdity 
of (A5). 



In analogy to (3.1) we define the coarse space, called further multiscale spectral coarse 
space, that is constructed using instead of {Cj}j=i by 

'tu ■= span{^j^^ \ j = 1, . . . ,nx and i = 1, . . . , Lj}, (5.5) 



where ipj are given by (2.4) with uiq. (•, •) replaced by rriQ^ (•, •) 



Remark 5.1. We note that for any subdomain Qj with Qj D d^l = 0, we have that (0, 'i-Uj) 
is an eigenpair of the generalized eigenvalue problem posed on Qj, where Iq^ denotes 

the constant 1-function on Qj. Thus, all multiscale partition of unity functions corre- 
sponding to subdomains flj with $7^ n dO, = are basis functions of our coarse space 'fn- 
This observation allows the interpretation of our method as a procedure that enriches a 
multiscale coarse space (given by the span of the multiscale partition of unity functions) 
by features that cannot be represented locally. These features are incorporated by those 
eigenfunctions corresponding to non-zero (but small) eigenvalues. 

6. Numerical Results 

6.1. General setting. In this section, we investigate the performance of the overlapping 
Schwarz domain decomposition method with coarse spaces discussed above when applied to 
some specific problems described in Section [4j We have implemented this preconditioner in 
C-|— |- using the finite element library deal.ii (cf. [2j). Our goal is to experimentally establish 
the robustness of additive Schwarz preconditioners with respect to contrast Kmax/'«min- 
The comparison is made using some coarse spaces known in the literature and the coarse 
spaces introduced in this paper. Namely, we consider the following coarse spaces. 

(1) Standard coarse space, denoted by := spanj^j | Cjldfij = for j = 1, . . . , Ux}, 
of standard partition of unity functions corresponding to interior coarse mesh 
nodes, which were introduced above (also, e.g [28|); 

(2) Multiscale coarse space, denoted by 'f'ff^ := spanj^j | Cjlan^ = for j 



1, . . . ,nx}, of functions that over each Qj are solutions of problem (5.1) and cor- 
respond to interior coarse mesh nodes (cf. |14j): 



(3) Spectral coarse space, defined in (3.1) as '■= span{^j(/3^ \ j = 1, . . . ,nx and i 
l,...,L,} 



(4) Multiscale spectral coarse space, defined in (5.5) (used only in the scalar elliptic 
case, see Section [s]). 

In all numerical examples (unless stated otherwise), the threshold for taking into account 
eigenpairs for the construction of the coarse space is chosen to be 0.5, i.e., I/ta = 0.5. 
Note that in the finite dimensional case (A5) is satisfied for any choice of the threshold 
1/t\. However, in practice one is generally interested in choosing 1/t\ not too large to 
avoid an unnecessarily large dimension of the coarse space. 



We consider Geometries 1 — 4 (see Figures 6.1 and 6.2), where for Geometries 1 — 3 k is 



equal to Kmin and Kmax in the white and black regions, respectively, and for Geometry 4 
K is given as shown in the logarithmic plot of Figure |6.2[b) . 
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0, 




^1 ^ 

Figure 6.1. Two sample geometries: Geometry 1 (left) and Geometry 2 
(right); the regions of low (white) and high (black) values of the coefficients. 
The mesh indicates the coarse triangulation. 



The geometries shown in Figure 6.1 differ by the number of connected subregions with 
high permeability. The goal for these two different distributions of the high contrast is to 

(1) test the robustness of the developed preconditioners with respect to the contrast and 

(2) show the benefits of the multiscale coarse space in the case of a large number of not 
connected, isolated, inclusions with high conductivity. 

In the abstract setting we replace the spaces 1^, ii), and Yh by their finite element 
counterparts. For this, we use the fine grid which is obtained from the coarse grid by 
subdividing the coarse grid elements into a number of finer elements. For Geometry 1 
and 2, we make an initial 8x8 mesh and introduce in each rectangular element an 8 x 8 
fine mesh, denoted by Th- Then, the spaces Yq, 'fo{^}j) and are finite element 

spaces corresponding to this mesh for a specific finite element, which needs to be chosen 
appropriately for the problem under investigation, e.g., Lagrange finite elements for the 
scalar elliptic problem in Galerkin formulation. In order not to overburden the notations, 
we have omitted the dependence of the spaces upon the fine-grid mesh size hoping that 
this will not lead to a confusion. 

Unfortunately, such choice of the spaces does not satisfy assumption (A4), since gener- 
ally the product (^jV of a standard partition of unity function and a finite element function 
does not belong to the finite element space Y. There are two ways to overcome this 
problem: (1) to project ^jV back to the finite element space using the form o(-, •) or (2) 
to use the finite element interpolant of in the finite element space. The projection 
is a local operation, but involves inverting some local stiffness matrices and could be un- 
necessarily expensive. The interpolation option, which was used in our computations, is 
straightforward to implement, but is not covered immediately by the abstract setting de- 
veloped above. However, a perturbation analysis could show that this is a viable practical 
approach whose rigorous study is a subject of our future research. 

6.2. Numerical Experiments for Geometry 1. In this subsection we use the standard 



coarse space Yff^ and the spectral coarse space defined in (3.1) 



6.2.1. Scalar Elliptic Problem - Galerkin Formulation (see Section 4-1)- Here, the finite 



element space is the space of Lagrange finite elements of degree 1. The right hand side / 



in (4.1 ) is chosen to compensate for the boundary condition of linear temperature drop in 
rr-direction, i.e., (f){x) = 1 — xi on di}. The dimension of the fine-grid space is 4225. In the 
PCG method we iterate to achieve a relative reduction of the preconditioned residual of 
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le— 6. In Tables 6.1 and 6.2 we present the results of two kinds of numerical experiments on 
the problem described in Subsection 4.1 for Geometry 1 with contrast Kmax/«^min increasing 



from le2 to le6. As partition of unity we use Lagrange finite element functions 

of degree 1 corresponding to the coarse mesh Th- 

In Table |6.1| we compare the number of PCG-iterations and the condition numbers for 
two preconditioners based on the standard coarse space (consisting only of the coarse 
Lagrange finite element functions) and the spectral coarse spaces Yh generated by our 
method, respectively. The standard coarse space has fixed dimension 49. The method 
performs well for low contrasts, but the condition number of the preconditioned systems 
as well as the number of iterations grow with increasing contrast. The spectral coarse 
space keeps the condition number independent of the contrast, which is in agreement with 
our theory. 

It seems that the results in the number of iterations for the space in the last row 
in Table 6.1 deviates from the general trend. We note that for all cases we run the 
PCG-method with the same stopping criterion, i.e., reduction of the initial preconditioned 
residual by a factor of le — 6. However, in this case the condition number of the precondi- 
tioned system is 1.77e5. Therefore, after reducing the initial preconditioned residual by a 
factor of le — 6 we may still be far away from the solution. Apparently, for larger condition 
numbers we may need many more iterations to compute the solution accurately. 





Standard coarse space 


Spectral coarse space i^n 




# iter. 


dim r^* 


cond. num. 


# iter. 


dim. fn 


cond. num. 


le2 


29 


49 


2.29el 


25 


76 


15.59 


le3 


50 


49 


1.88e2 


21 


145 


11.51 


le4 


55 


49 


1.79e3 


18 


162 


6.20 


le5 


67 


49 


1.78e4 


18 


162 


6.18 


le6 


66 


49 


1.77e5 


19 


162 


6.19 



Table 6.1. Elliptic Problem of Second Order: Numerical results for stan- 
dard and spectral coarse spaces 



In Table 16.21 we show the number of PCG-iterations and condition numbers for two 
preconditioners based on spectral coarse spaces. In columns 2 — 4 we report the results 
for a coarse space of fixed dimension 162 and the threshold for which this is achieved. In 
columns 5 — 7 we present the results for a fixed threshold 1/t\ = 0.5. We note that the 
difference in the performance is only for values of the contrast below le4. 



6.2.2. Scalar Elliptic Problem - Mixed Formulation (see Section Here, the finite 

element space is the Raviart-Thomas space of degree (RTO) for the velocity and piecewise 
constants for the pressure on the same rectangular fine mesh as above. The right hand side 



/ in (4.6) is chosen to compensate for the boundary condition of unit flow in x-direction. 



i.e., u ■ n = ei ■ n on dil., where ei is the first Cartesian unit vector. The (divergence free) 



coarse velocity space is constructed as outlined in Remark 4.1 We first construct a basis of 
the spectral coarse space corresponding to the stream function space. The corresponding 
coarse velocity space is then given by the span of the curl of these basis functions. Note, 
that the stream function space corresponding to RTO is given by the space of Lagrange 
polynomials of degree 1 (see [THl Section 4.4]). 
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Spectral coarse space of dim "Vn = 162 


Spectral coarse space, 1/tx = 0.5 


'^max/ ^min 


# iter. 


cond. 


I/ta 


# iter 


cond. 


dim fn 


le2 


18 


7.45 


1.39 


25 


15.6 


76 


le3 


17 


5.99 


0.92 


21 


11.5 


145 


le4 


18 


6.20 


0.50 


18 


6.20 


162 


le5 


18 


6.18 


0.50 


18 


6.18 


162 


le6 


19 


6.19 


0.50 


19 


6.19 


162 



Table 6.2. Elliptic Problem of Second Order: Numerical results for spec- 
tral coarse spaces of fixed dimension and fixed threshold 1/t\ = 0.5. 





Standard coarse space 


Spectral coarse space i^n 


'^max/^min 


# iter. 


dim rj* 


cond. num. 


# iter. 


dim -Th 


cond. num. 


le2 


32 


49 


2.87el 


23 


86 


13.87 


leS 


50 


49 


2.26e2 


24 


129 


18.38 


le4 


63 


49 


2.19e3 


17 


162 


6.57 


le5 


80 


49 


2.18e4 


18 


162 


6.65 


le6 


87 


49 


2.13e5 


19 


162 


6.68 



Table 6.3. Scalar elliptic equation in mixed formulation: numerical re- 
sults for Standard coarse space and spectral coarse spaces obtained 
for a fixed threshold 1/tx = 0.5. 



As partition of unity we could simply use the bilinear Lagrange basis functions 

corresponding to the coarse mesh Th- Nevertheless, for consistency with the Brinkman 



case (see Section 6.2.3), where we have higher regularity requirements, we choose the ^^'s 
to be piecewise polynomials of degree 3, such that all first derivatives and the lowest mixed 
derivatives are continuous and Cji^i) = for i,j = 1, . . . ,nx- 

In Table 6.3 we present the numerical results for this problem and Geometry 1 (see 
Figure 6.1). The dimension of the fine space is 12416. In columns 2 — 4 we report the 



number of iterations, the size of the standard coarse space, and the condition number of 
the preconditioned system. Here, the standard coarse (velocity) space is given by the span 
of the curl of the partition of unity functions corresponding to interior coarse mesh nodes. 
Columns 5 — 7 contain the number of iterations, the dimension of the coarse space, as well 
as the condition number of the preconditioned system. It is clear that for the standard 
coarse space of dimension 49 the condition number grows with increasing the contrast 
and so does the number of iterations. However, when the coarse space includes all coarse 
eigenfunctions below the threshold 1/tx = 0.5, the preconditioner shows convergence rates 
and condition numbers independent of the contrast. 



6.2.3. Brinkman Problem (see Section 4-4^- Next, we present the numerical experiments 
for the Brinkman problem (4.13), where the right hand side / is chosen to compensate for 
the boundary condition of unit flow in x-direction, i.e., u = ei on d^l. The viscosity n is 



chosen to be 0.01 and k varies depending on the contrast (see Table 6.4). 
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Standard coarse space "V^ 


spectral coarse space 'Vu 


^max/'i^min 


# iter. 


dim. 


cond. num. 


# iter. 


dim. 


cond. num. 


le2 


27 


49 


2.13el 


25 


60 


14.69 


le3 


39 


49 


4.22e2 


28 


75 


21.73 


le4 


70 


49 


2.25e3 


29 


106 


21.83 


le5 


91 


49 


1.51e4 


24 


153 


13.08 


le6 


113 


49 


1.24e5 


22 


164 


13.82 



Table 6.4. Numerical results for Brinkman's equation using standard 
coarse space and spectral coarse spaces 



We discretize this problem with an iJ((iiu)-conforming Discontinuous Galerkin dis- 
cretization (cf. [m [32]) using Raviart- Thomas finite elements of degree 1 (RTl). We 
again employ a 64 x 64 fine grid. It is well-known (see [131, Section 4.4]) that in two spatial 
dimensions the stream function space corresponding to the RTl space is given by Lagrange 
biquadratic finite elements. For a generalization to three spatial dimensions one has to 
utilize Nedelec elements of appropriate degree. As above, we use an 8 x 8 coarse mesh. 



We choose \^^jY\=\ described in Section 6.2.2, which satisfies all regularity constraints. 

In Table |6.4| we give the number of iterations, the dimension of the coarse space in 
the additive Schwarz preconditioner, as well as the estimated condition number of the 
preconditioned system. The dimension of the fine space is 49408. As for the scalar elliptic 
case in mixed formulation, the coarse (divergence free) velocity space is constructed as 



outlined in Remark 4.1 In columns 2 — 5 we present the results for the case of the 



standard coarse space of dimension 49, which as in Section 6.2.2 is given by the span of 
the curl of the partition of unity functions corresponding to interior coarse mesh nodes. 

We observe, that the increase in the contrast leads to an increase in the condition 
number and subsequently to an increase of the number of iterations. Further, in columns 
5 — 7 we report the number of iterations, the dimension of the coarse space in the additive 
Schwarz preconditioner, as well as the estimated condition number of the preconditioned 
system for the spectral coarse space obtained by a fixed threshold \It\ = 0.5. For the 
Brinkman problem the performance of the preconditioner is also robust. We should note 
however, that Brinkman's equation is much more difficult to solve due to the fact that the 
overall system of linear equations is a saddle point problem. 



6.3. Numerical experiments for Geometry 2 in Figure 6.1 These numerical ex- 



periments are aimed to compare the performance of the iterative method applied to the 



second order elliptic problem in Galerkin formulation (see Section 4.1) for a permeability 



given in Geometry 2 (see Figure 6.1). The goal here is to demonstrate the coarse space 
dimension reduction when using multiscale partition of unity functions instead of standard 
ones. The dimension of the fine-grid space is 4225. 

In Table |6.5| we present the results when the preconditioner is based on the spectral 
coarse space (columns 2 — 4) , the multiscale spectral coarse space Yh (columns 5 — 7) , 
and the multiscale coarse space Yff''^ (columns 8 — 10). Comparing the data for the 
spaces 'Vh and shows that the number of PCG-iterations and the estimated condition 
number of the preconditioned system are robust with respect to the contrast Kmax/^min- 
We can also observe that when using the spectral coarse space the dimension of the 
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'^min 


^ iter. 


dim Th 


cond. # 


# iter. 


dim 


cond. # 


# iter 


dim 


cond # 


le2 


22 


163 


12.15 


21 


44 


10.81 


19 


49 


8.70 


le3 


18 


612 


8.42 


20 


60 


9.86 


35 


49 


5.97el 


le4 


15 


838 


4.92 


22 


60 


10.90 


44 


49 


5.63e2 


le5 


16 


838 


4.92 


22 


60 


11.01 


53 


49 


5.59e3 


le6 


17 


838 


4.92 


22 


60 


11.01 


66 


49 


5.59e4 



Table 6.5. Scalar elliptic - Galerkin formulation: results for spectral 
coarse space i^H , multiscale spectral coarse space Yh , and multiscale coarse 
space r^". 



coarse space increases as the contrast increases Kmax/'^min, which is in agreement with the 



analysis of Section 4.1 The decrease in the estimated condition number when going from 
/^max/Kmin = le2 to Kmax/'^min = le3 and further to Kmax/'^min = le4 Can be explained 
by the fact that for higher contrasts more eigenvalues are below the prescribed threshold, 
yielding a higher dimensional coarse space and a lower condition number. However, it is 
important to note that the dimension of the coarse space reaches a maximum for Kmax/«^E 



above a certain threshold. As we can see for /^max/^min In the range l64, . . . , l66, the 



dimension of the coarse space stays the same. By the analysis in Section 4.1 we know that 
there is only a finite number of asymptotically small (with the contrast tending to infinity) 
generalized eigenvalues. The reported data can be seen as evidence that for this specific 
configuration we have reached this asymptotic regime for Kmax/^min = le4. 

In columns 8 — 10 we present the numerical results of the algorithm when the precon- 
ditioner is based on the multiscale coarse space '^ff^, which consist of one basis function 
per interior coarse node. As we can see from the data, the number of PCG-iterations as 
well as the condition number of the preconditioned system grow steadily with the growth 
of the contrast. ^ 

The important point to observe when using the multiscale spectral coarse space 'Vh 



(see columns 5 — 7 of Table 6.5) is that its dimension is drastically reduced compared 
to the spectral coarse space Yh- In our specific example, the largest dimension of the 
multiscale spectral coarse space 'fn, which is constructed using the multiscale partition of 
unity {^j}"^]^, is 60, compared to the dimension 838 of the spectral coarse space Yh, which 
is based on the standard partition of unity {Cj}j=i- One is generally interested in keeping 
the dimension of the coarse space as small as possible, especially when the problem is 
solved multiple times. The data is a confirmation of our reasoning in Section [5] 



6.4. Numerical experiments for Geometries 3 and 4. In Table 6.6 we present the 



numerical results for the scalar elliptic equation of second order in Galerkin formulation 



from Section 4.1 for highly heterogeneous permeability distributions shown in Figure |6^ 
Geometry 3 represents a rather challenging example: the permeability field is highly het- 
erogeneous with more than 4000 small and about 100 large randomly distributed inclusions. 
We consider this a challenging test for the robustness of the iterative method by perform- 
ing a relatively small number of iterations using a coarse space of low dimension. Here, 
we have used a 16 x 16 coarse mesh and subdivided each coarse cell into 16 x 16 subcells 
to obtain a 256 x 256 fine mesh. The preconditioner is based on the multi-scale spectral 
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Figure 6.2. Two geometries with a 256 x 256 fine mesh and a 16 x 16 coarse mesh. 



coarse space Yh- The dimension of the fine space is 66049 while dimension of the coarse 
space is at most 293. As we can see, the condition number of the preconditioned system is 
robust with respect to the contrast and the dimension of the coarse space is quite small. 

Geometry 4 (see Figure 6.2 ) represents a more challenging problem in that it is no longer 
a binary medium, i.e., k assumes many and not just two extreme values. The geometry 
is generated by setting k in a fine mesh cell to IQT'''^'^'^, where rand denotes a uniformly 
distributed random number and 7 = 2, . . . , 6. This produces a random field 10'^^^^^ where 
r]{x) is a realization of a spatially uncorrelated random field. This yields a "background" 



on top of which we put randomly generated inclusions similar to Geometry 3. In Table 6.6 



(columns 5 — 7) we report the numerical results using the preconditioner based on the 
multiscale spectral space ^h- As we can see, the number of PCG iterations and the 
condition number of the preconditioned system are robust with respect to increases in the 
contrast. It is furthermore important to note that, even for this random case, the dimension 
of the coarse space stays reasonably small (at most 397) compared to the dimension of the 
fine space, i.e., 66049. This exemplifies the robustness and applicability of the numerical 
method developed above. 



7. Conclusions 

The theory developed above introduces a method for constructing stable decomposi- 
tions with respect to symmetric positive definite operators. The robustness with respect 
to problem and mesh parameters is proved under rather general assumptions. We have 
furthermore applied this abstract framework to several important cases, i.e., the scalar 
elliptic equation in Galerkin and mixed formulation, Stokes' equations, and Brinkman's 
equations. For the scalar elliptic equation in Galerkin formulation, we have additionally 
presented a strategy of reducing the dimension of the coarse space in the stable decompo- 
sition. To verify our analytical results, we have performed several numerical experiments, 
which are in coherence with our theory and show the usefulness of the method. 
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Geometry 3 


Geometry 4 


'^max/ ^min 


# iter. 


dim >^ 


cond. # 


# iter. 


dim -fu 


cond. # 


le2 


22 


209 


11.2 


19 


217 


8.47 


le3 


24 


259 


15.2 


20 


221 


9.48 


le4 


23 


275 


11.2 


22 


244 


11.5 


le5 


24 


277 


11.2 


25 


317 


16.3 


le6 


27 


293 


11.7 


23 


397 


11.7 



Table 6.6. Scalar elliptic problem - Galerkin formulation: Numerical re- 



sults for permeability fields shown in Figure 6.2 using multiscale spectral 
coarse spaces 'fn- 
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